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SUMMARY 


The heat transfer due to catalytic recombination of a partially 
dissociated, diatomic gas along the surfaces of two-dimensional and 
axisymmetnc bodies with finite catalytic efficiencies is studied ana- 
yt ic ally. An integral method is employed resulting in simple yet 

re a ively complete solutions for the particular configurations 
considered. 


A closed form solution is derived which enables one to calculate 
atom mass-fraction distribution, therefore catalytic heat transfer dis- 

If'- ^ he surface of a flat plate in frozen compressible flov, 
with and without transpiration. 

,, _ Numerical calculations are made to determine the atom mass- fraction 
distribution along an aicisymmetric conical body with spherical nose in 
frozen hypersonic compressible flow. A simple solution based on a local 
similarity concept is found to be in good agreement with these numerical 
calculations. The conditions are given for which the local similarity 
solution is expected to be satisfactory. 

fl * „ h + h nf' 1 i!! lt S i0 f the practical application of the analysis to the 
light of the blunt bodies m the atmosphere are discussed. The use of 
boundary-layer theory and the assumption of frozen flow restrict appli- 
ca ion of the analysis to altitudes between about 150,000 and 250,000 


INTRODUCTION 


. Surface reaction, long a subject of interest in the chemical 
industry has become important in aerodynamics in connection with hieh 
a 1 ude hypersonic flight. When the flight conditions are such that 
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_ oir -n^t the body is almost chemically frozen , 

the flow of dissociated aip influenced by the amount of cata- 

heat transfer to the body is f 011 ® J/" lYti ; wall . the chemical 

lytl0 ; r r "S transf^ at aU an,, 

energy of the flow ^es c val] _ s the heat transfer is approxi- 

conversely, with highly ^ across ar. equilibrium boundary 

mately the same as that wh , ,n <?qr dated the heat transfer 

layer, ttus, when the air is =°" P 'f ^the ^ if ffom about 25 percent 
may vary with the catalytic activity ottne w oi 

to about 100 percent of equilibrium heat ran:. 

, nr ,„it Pf ] n surface chemical reaction 

Boundary-layer theory was first app ' people, such 

problem by Chambr4 and f ™os <ref. l).Ji £ 

through^ 8^ lave* studied theVoblem for varlois appllcatlone using dif- 
fered approaches. The specific solutions obtained are for 

1 . ®e incompressible flow over a flat plate upon which an arbi- 
trary order chemical reactioni OOCUIB (refs. 1 through 5 ). 

2. The compressible flow over a flat plate upon which first-order 
chemical reaction takes place (ref. ). 

o The first-order chemical reaction at the stagnation point m a 
3 - dissociated hypersonic flow without mass transfer (ref. o) , and 
with mass transfer (ref, 7)* 

T _ +his Ba p e r an integral method will be used to obtain catalytic 
In this paper, an n s i e in frozen compressible 

heat-transfer predictions for ( ) axisymmetric conical body 

s rss »« -**“ 

transpiration. 

The effect of mass transfer is Includes in the present 

rsrctl - surface 

participates in the chemical reaction. 


A 

)| 
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SYMBOLS 


Bl ,B 2 constants for equation (51) given In table, page 16 

/ Pw.o^oY ' 2 
L \Pe,o^e,o/ 


OThe reaction rate at the wall is an arbitrary function of the 
surface concentration of the reactant. 
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C 

c 

C P 

D 

F 

Fi 

F 2 

G 

G 

h 

Ah° 

J V 

K x 

Ko. 

K 2 

k 3 

k 

k w 

L 

M 

m 

P 

Pr 


parameter defined, by equation (9) 
mass fraction 

frozen specific heat at constant pressure 
coefficient of binary diffusion 
Blasius stream function used in equation (63) 
momentum integral defined by equation (39) 
concentration integral defined "by equation ( 4l) 
reaction parameter defined by equation (49) 

G 

n/Po 

frozen total enthalpy defined by equation (6) 
heat of recombination 
mass flux of atoms to the wall 
parameter defined by equation (53) 

Ja 

parameter defined by equation (54) 

parameter defined by equation (55) 

frozen thermal conductivity of gas 

specific rate constant for surface recombination 

reference length (see fig. l) 

molecular weight 

free- stream Mach number 

c i 

c i,e 

pressure 
Prandtl number 





q 

lc 

^d 

R 

r 

Re 

s 

Sc 

T 

t 

U 

u 

V 

v 

X 

x 

y 

z 

z 

3 

r 


total heat-transfer rate toward the wall per unit area 

convective heat-transfer rate toward the wall per unit area 

heat transfer rate per unit area hy atoir recombination at the 
surface 

universal gas constant 

distance from the axis of symmetry to a point on the surface 
(see fig. l) 

. -i PTJqqL 

Reynolds number , — — 

st Teamwise independent variable defined by equation (T) 

|JL 

Schmidt number, 
absolute temperature 

independent variable defined by equation (8) 

_u = Si 

U e Lit 

x component of velocity 
Os 

y component of velocity 

2 L 

L 

streamwise distance along the surface 
distance normal to surface 

mass rate of collision of atoms with tie wall 

variable defined by equation (65) 

due 

dx 

transpiration parameter defined by equation (38) 



surface catalytic efficiency 
boundary-layer thickness 

modified boundary-layer thickness defined by equation (36) 

JO for t vo- dimensional bodies 
(1 for axi symmetric bodies 

jrj; defined by equation (36) 
nose angle 

parameter defined by equation (37) 

dynamic viscosity 

kinematic viscosity 

dummy variable 

density 

m v _ c i,v 
m w^o Cl,V,0 

shape factor defined by equation (50) 
shock layer thickness 

o $ 

equilibrium mole fractions 

stream function defined by equations (l4) and (15) 
parameter defined by equation (8l) 

Subscripts 


equilibrium 

edge of boundary layer 
ith species 



6 


m mixture 
w wall 

0 leading edge of flat plate or stagnation point of sphere 

1 atoms 

2 molecules 

oo free stream 


Superscript 

differentiation with respect to s except for equation (63) 


A 

h 

2 
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SIMPLIFICATION OF BOUNDARY- LAY! R EQUATIONS 


The classical laminar boundary- layer theory is used. in the present 
paper. The fluid is considered to be a partially dissociated diatomic 
gas of constant Schmidt number and it is assumed to be frozen in is 
state throughout the entire flow field. The actual range of lgh^ 
ditions in the atmosphere for which these simplifications aie valid 
be discussed later. 

The equation of state for a partially dissociated diatomic gas is: 


P ' | 


(1) 


where the term (l+c x ) is known as the compressibility factor. For wo 
dimensional and axisymmetric bodies the bounlary-layer conserva ion 
equations are: 


Mass 


Hpur 6 ) d(pvr e ) _ Q ( 2 ) 

■x oy 


dU 

pu — + pv 

OX 


Ju 


du 

= — l ^ — 


jy dy 




Momentum 


_ SB 

dx 


(3) 
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A 

), 

2 
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Frozen total energy 


pu ^ + pv ^ = 
dx dy 




0) 


Species 




it Aj_ ijcA 

d V Sc dy J 


(5) 


where 


h = 



+ 


2 


(6) 


Considerable simplification of the conservation equations results when 
a new set of independent variables is defined as 


s 


Cu f >r 2€ dx 


(7) 



Following references 9 and 10 we assume 

P _ f °w,o^tf : o V 

V^ejO^ejO/ 

and 

Pe^e _ Pe 
Pe,oP-e,o Po 

Now C is a function of x alone . It was shown in references 9 and 10 
that these assumptions (eqs. (10) and (ll)) affect the solution for heat 
transfer to an infinitely catalytic wall or a totally noncatalytic wall 
by less than 10 percent. It might be expected that discrepancies for 
intermediate values of catalytic activity will be as small. 


( 10 ) 

(ID 
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Transformation of the conservation equations to the set of 
independent variables s and t is accomplished by means of the following 

relations: 



Cu e r 2€ 


M2 + at a(__) 

Ss 8x }t 


( 12 ) 


A stream function is 
(2) as 


M2 , u r* 

Sy e Pe,o 
defined to satisfy the nass 


( 13 ) 


conservation equation 


pur G 


pvr 6 


= P 


= -P 


St 

(14) 

Sy 


St 

(15) 

Sx 



A 


Transformed to the s— t plane , equations (id) and (15) become 



u. = St 
"■J- e St 

( 16 ) 

pv = 

- £ 0 ^ fouer* & + 4* $) 

r e \ ^ 5 s Pt/ 

( 17 ) 

Nev dependent variables 

are defined as 



v-g 

5t 

( 18 ) 


V = - §S 

( 19 ) 


5s 


Now with the use of equations ( 7 ) through (3-9), the conservation equa- 
tions (2) through (5) can be transformed to somewhat resemble the 
familiar incompressible boundary- layer equations as follows: 


Mass 


Sy + Sy = 0 

8s St 


( 20 ) 


OJ LT\ 



< -3- OJ ITN 


2G 


9 


Momentum 




os 


St 


c)t 2 u e ds P 


Frozen total energy 

U _j_ y Sh _ y o S h 

Species 


s £ 


St Pr St 2 


+ v, 


> , e u e ( 1 


1 ) S 2 At 


Pr J St 2 V 2 


Sc * Sc.. V, 

U — i + V — ± = -±U 

os dt Sc ^t 2 


o,e ^ 2c i 


( 21 ) 


( 22 ) 


(23) 


The last. term in equation (21) is zero for a flat plate without 
pressure gradient. Also , according to reference 9, it may be neglected 
in a highly cooled hypersonic boundary layer over a blunt body. Equa- 
tion (21) therefore becomes for a flat plate or a highly cooled hyper- 
sonic blunt body 


u2U + v2 !J 

OS St 


'o , e 


Sfu 

ot 2 


(24) 


Equation (23) may be written, for the 


diatomic gas under consideration 


u 2£i + v 

os 


Sc x 

St 


v o , e S 2 c 1 
Sc St 2 


(25) 


C 2 - 1 - Cjl ( 26) 

Note that as a consequence of the assumption of equation (10) the 
momentum and species conservation equations are no longer coupled to the 
energy equationj therefore, equation (22) need not be treated here since 
its solution may be found elsewhere. Equations (20), (24), and (25) are 
in the form of the familiar incompressible flat plate boundary-layer 
equations and are now amenable to solutions. Now it is seen that the 
pertinent original boundary-layer equations (2), (3), and (5) are simpli- 
fied to equations (20), (24), and (25). The boundary conditions on these 
equations are as follows: 
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at t = 0: 


U = 0 , v = 


Po,e Cu er e 


d Cl _ Po.e Sc 


Sc VvjCi.w Sc T 

— 2 — + — Jx w 

Po,e P-o,e^ u e r "' 


U = 1 , 


; i = c i,e 


sTsim-* Bounded , ^ (§J 


Bounded. 


“\rf ary r^iHhnasffS o^oms tfthe ^ aid T 
known. The term J i? w * chemical kinetics of the surface 

cbdained ^ “ * hall include in the iubsequent section a brief 

SscusSh S "he surface chemical reaction ki reties before going into 

the actual solution of the equations. 

CHEMICAL KINETICS OF THE CATALYTIC SURFACE REACTION 


More comprehensive treatments of the sub,; 
may be found in references 11 and 12. 

It is generally accepted that the catalyt 
dated air on surfaces at temperatures lower t 
above 10- 4 atm is a first-order process procee 
ligibly small equilibrium atom concentration, 
to such a reaction* 

The microscopic reaction rate, which is 
collision with the wall, may be obtained as f . 
diatomic gas: Let 7X 2 Tie the probabi 1 Y ° 
striking the wall. The equilibrium mole frad 
is unity at the surface temperatures and pres: 
The factor y is commonly referred to as the 
The mass rate of recombination may then be wx 


ect about to be discussed 

ic recombination of disso- 
,han 2000° K and at pressures 
;djig only toward the neg— 
This discussion is limited 


i function of the rate of 
allows for a dissociated 
f an atom recombining upon 
tion of the molecules, X 2 , 
sures of interest here. 

"catalytic efficiency." 
itten 



< OJ LT\ 
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where Z is the mass rate of collision of atoms with the wall and is 
for a Maxwellian distribution, 


then 



(30) 


(3D 


In engineering work, including the present analysis, the reaction 
rate is often written 


J i,w = kwPw^i.w = 

The specific catalytic rate constant 
efficiency by the expression 


2k. 


PMx 


c i^w 


, k w , is related to the 


(32) 


catalytic 




RT- 


w 


V W 


(33) 


TIig magnitude and temperature dependency of 1 ky ape not well 


known for most surfaces; however, metallic and metallic oxide surfaces 
have much higher efficiencies than nonmetallic surfaces. For instance, 
at room temperature, k w for a metallic oxide is about 10 ft/sec, whereas 
ky Pyrex is about 0.1 ft/sec (see ref. 6). There are also possi- 

poisoning" surfaces to reduce their catalytic efficiency. 
Reference 13 reported a substantial reduction of the catalytic efficiency 
of metallic surfaces upon application of iodine coatings. 


REDUCTION OF GENERAL EQUATIONS TO AN ORDINARY DIFFERENTIAL EQUATION 


In this section, the partial differential equations (20) , (24), and 
( 25 ) with the boundary conditions of ( 27 ) will be reduced to an ordinary 

equation by the use of an integral method described in ref- 
erence 14. 

The diffusion boundary layer is assumed to have the thickness of the 
momentum boundary layer, an assumption which is believed to be valid for 
Sc in the order of 1 according to reference l4. Equations (24) and 
( 25 ) are integrated across the boundary layer with the aid of equation 
(20) . The resulting integro-differential equations are 
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where 


AF 


t 

1 




+ r 


AF 2 ' + ^ A' 






A = 


&t 2 

v e ,o 


VySt 

Ve ; o 


Fi = 


U(l-U)dq 


m 


Cl 

c i,e 


F 2 = / U(l-m)dr) 

J o 

and the prime denotes differentiation with reject to s. The 
conditions become, with the aid of equation (,3 )> 

at T) = 0: ' 

U = 0 


V = v w 



r Sc in + 


pk v Sc & t 


m 


> 


j 


(3l) 

(35) 

(36) 

(37) 

(38) 

(39) 

(10) 

( 11 ) 

boundary 

( 12 ) 


<4-4- OJ ITS 



CVJ LT\ 


13 


at t] = 1 


U = m = l 

^ = ^ = 0 
dT i 


03 ) 


mst Sr^ he e ^ atl ° ns (34) and (35), the form of U and m profiles 
must first be assumed. Fourth and fifth degree polynomials in / are 

chosen to represent the profiles for U and m, respectively. 

Thus: 


U — • ) 3^T]' 
11=0 


n 


om 


m = 


) t»nTi n 


(4-5) 


n=o 


The coefficients a n and b n are obtained by satisfying a sufficient 
s^liv/f*/ COTdlti °" s « 1 = 0 an/, = 1. 

lef/tn l , ^ ountol T conditions. One coefficient in equation (/) is 
left to be determined as the solution of equation (35); therefore only 
five boundary conditions are needed. These boundary oindit(ons are 
those from equations (te) and (k 3 ) , and those obtained by saSsmL 

“ “ d \T ° 1-0- The particu/r/hoicTof 

profiies, equations (44) and (45), is justified a posteriori by the 

as wilZb/show/lat er . present -^tions with available exact solutions 

“fth“ pre- 


A' = 


( OU/orj ) 


W 


+ r 


Fx/2 


= constant 


(46) 


Therefore: 


A = A 1 s 


( 47 ) 


l4 


The catalytic bounty conaiticn descHW in («) W he -men, 
with the aid of equations (l) , (T), (9). ^ < f) > 


/-v \ UlW 

= r Sc mw + G Sc 0 T+mwCi,e 


(48) 


where the reaction parameter, G, 


IS 


G = 


(l+ei.e)Te,o 

T w 




m 


and the shape factor , 0 , is 


$ 


/ 


CLu e 2 r 2e 



(50) 


Upon substitution of the ( 4l)^ becomes 1 

boundary conditions, the concentration integral, v ), 


F 2 — Bi(l-%) — 



(51) 


where the constants 
cnlutiion of 


B x and B 2 depend on 

'Lt.inn (Ik) incorporated 


Equation (35)> vith the 
now he written as 


K 3 - K 2 m w - kJ 

/ 2 ^ 
^ + A'B 2 Scy 

| q + 2$ 1 s j 

m w 

| 1+mwC i,e 

2 I 

K 2 + Ki<& > 2 

(l+mwC;.^e/ J 

l S 


(52) 


where 


Kl - ¥ 1 G So 
Ka - 1 + Fa r S: 
Ka = 1 " x% 1 


(53) 

(5*0 


( 55 ) 


<; cvj ir\ 


^Jl l\J 
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It is now seen that the pertinent 
reduced to the single equation (52) . 


conservation equations are 


SOLUTION FOR THE SURFACE DISTRIBUTION OF ATOM CONCENTRATION 

Condition for Application of the Local 
Similarity Concept 


A 

4 


Before actually solving the equation (52), let us examine the pos- 

S1 1 * ^ 0 obtaining the surface distribution of atom concentrations by 
use of the local similarity concept. J 

, . f } OW p f st a wal1 with a finite catalytic activity, a concentra- 

lon boundary layer develops through which atoms diffuse to recombine on 
e surface. The resulting steady state atom concentration along the 
wa is determined by the relative local rates of atom diffusion to the 
ur ace and catalytic recombination at the surface. The diffusion rate 
oi the atoms, for a given concentration potential, is determined by the 
aerodynamics of the flow, characterized here by the momentum and diffu- 
sion boundary- layer equations ( 24 ) and (25). The consumption rate of 

chel S S " al U ° n the ° ther h£md ' iS control l e <l by the heterogeneous 
chemical reaction kinetics discussed earlier. In general, affine solu- 

10ns of equation (25) can not be obtained to satisfy the catalytic 
boundary condition of equation (27). This can be seen, for instance 
10 ™ e f act that the catalytic boundary condition ( 18 ) for the inte- 
grated form of the diffusion equation, (35), depends upon the variable 
X explicitly through According to the present method of solution 

characterized by the equation (35) and its boundary conditions, one would 
have a similarity solution which would specify that m w be independent 
_X for any body whose geometry and flow characteristics are such that 
is independent of X. On the other hand, when <t varies slowly with 
, one expects a local similarity solution, explained below, to be valid. 

Here a local similarity solution means that equation (35) is solved 
localiy with the boundary condition ( 18 ) applied as though it were not a 
unction of X, however, with 0 being calculated at the point in ques- 
tion. An affine solution can in this way be obtained locally for each 
position along the body. 

In practice, one could make use of already existing affine solutions 
for a thermal boundary layer over an isothermal flat plate (ref. 15 for 
instance) . This problem is analogous to the present one when the atom 
concentration at the nail is independent of X. The affine solution 
yieids a relationship between (dm/c-T)) w and m w for a locality which is 
then solved simultaneously with equation ( 48 ) for m w (x) . 
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The hypersonic pressure distribution and equilibrium heat-transfer 
already ^ 

t ualitatively 1 checked Stjlout performing any calculations. The simple 
qualitati y _ rpference 15 for instance, describes quite accu- 

“"illbrt™ heat transfer to t-dimenslonal 
5££r obt2n 

( 56 ) 


lw,E 

c lw,o,E 


Pe 1 
Po ® 


=r £S.r;:: 

are found to be “^^distribution of atom concentration 

‘Sng'Se wall' therefore the catalytic heat transfer, may be obtained 
simply by using the local similarity concept. 

-CtiSs 

The universal constants involved ^Tre"" 

the three transpiration rates, r - 0, l.u, » 
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r 

X ' 

Bi 

B2(f° r Sc = 0.72) 

0 

34.05 

0.2429 

0.06270 

1.0 ] 

45.64 ] 

.2323 

.06457 

2.0 

58.51 

.2244 

.06679 


Flat Plate With Similarity Transpiration 

The body studied in this section is shown in figure l(a). For two- 
dimensional flow without pressure gradient, the streamwise variable 

becomes 


s = CLueX 


(57) 


<3-4* OJ 




< -=* cm ir\ 


G 


$ =Vx 


Equation ( 52 ) becomes 


5£w _ k 3 - **» - ( 2 + 1^—7 

dx : 1 ( 

o V K x x 2 

2 K 2 + ± X 

L ( 1 +myC 1 ^ e ) 2 J 

inlls Ug ltlTnZ:Vej “? 8mted readily ’ md these '“<* “ill »* 

1. When the surface is infinitely catalytic. 

2. When the surface is totally noncatalytic. 

3 * Wh “- the com P ress ^ility factor at the wall, (l+m^ 0 ) is 
sufficiently near 1 . ’ K m u<-i,e) , is 


For the first case it is obvious that 


m w - 0 for all X. When the wall 


is noncatalytic, the similarity solution^!' equation ( 59 ) j 


m ir = — 

w Kp 


near ^ 5 ^ ™ U ' , is sufficiently 

’ ^ (59) simplifies to a first-order linear equation 


d^w V A’B 2 ScJ 

dX TT— 

L 2(K 2 +KxX2)x 


1 

K x X 2 


m v = 


2(K 2 +K x X 2)x 


Swonrtiat 01 ^ 10 ^'?^ :ri“f satisfies the — 


1 - fi + |i x iVv 1 + 


2 

A*B 2 Sc 


1 + — 2 _f\ £1 x i 
VB a Scy k 2 
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A limiting process , letting X -* 0 , 
result since at X = 0 there is no 
inert wall solution, equation (60), 


shows that mw,o = K3/K2, an expected 
time for reaction. Therefore the 
must be v^lid. 


More exact solutions of the limiting g^al^compare^hese solu- 

( &o\ -Pound elsewhere and we snaxx corny cii 

tions (60) an ^^ a ^f s ^ uti on for noncataly.ic wall with similarity 
tions here. The exact soiuuu reference 7 as 

transpiration can he derived from equation (1» of referen- 


% = 


1 - F-ySc 


,n 

fell] 

Jo 

|_F"(°)_ 


ai 


( 63 ) 


The transpiration rates are related by. 


= - 2 


v w I p-w^gX _ _ 2f 
Ue n/ hw ./v 


(64) 


_ _o.5 (separation occurs at F w 

_ . 1 • -1- ^ nnnof 1 OP 


- 2385 ) and Sc = 0. (2, the 

For i' = -u.p iq ui thin 6 percent of the 

rarSuticn^rvSuHf ^"integral in aquation (63) was obtained 

from reference l6. 

»ow ton the case where the cOBpressibilLty facton at ^ is 

sufficiently near 1, T - 0, and references 1, 4, and 5, 

££ compressibility and ^s^raUon were not tahen into account, 
terras of the parameter z used in reference 1, 


In 


z = 


k w L Sc p e 

0.339^e 


Re 


1 

2 Sc 


X 1 

' :i x 2 


(65) 


the present solution for Sc - 0. f2 is 


_ 1 - (1+0.458z)~ 
m ' v ' = 1.053z ” 


•30 


( 66 ) 


, . Q+ . n n t(jc\ is compared with other known solutions in 

This solution, equation (6 >> ement is shown between the present 

IZZolSZlt of°references 1 and 

sta£“ue to Z STS aXffe^nt orSer polynomial approximation 
in the velocity profile. 


VJ1 4^ > 
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A 

4 

2 

5 


The complete solution of equation ( 59 ) is obtained numerically by 
the use of an IBM ^digital computer. The numerical solutions were 
o tamed for Sc - 0.72 and several combinations of K x and c, p at 

transpiration rates of T = 0 10 and 2 D A o +11 ,n r 1 ' e . 

+ Q „ ’ v > ana A study of these numerical 

resuits revealed that for each value of r, the solutions were congruent 

(6 ?a Vh6n 0 = m w/ m v,o Was P lotted versus § 

lX /K 2 (l+c 1? e) CT . Therefore, the function K x X 1/s /Ks in equation (62) 

TllslTlT ed b ( t KlX and the follovlng^eMe-pirLa? ‘ ’ 

rm solution was obtained for all compressibility factors. 


a = 


m w 


m. 


1 - 

r, ^x 2 

1 + 1 1 

1 - 1 + - - 

\ A'B 2 Sc 


K 2 (l+Ci^ e ) a J 



1 + 


i 

K x X 2 


( 67 ) 


A'B 2 Sc J K 2 (l+c l7 e) a 


Figure 3 shows the close correlation between equation (67) and the 

results of the numerical solutions of equation (59) for the three trans- 
piration rates „ 

It 1S seen in equation (67) that the general solution depends on 
the two parameters, K a X 1/2 /K 2 (l+c^ e ) a and 2/A’B 2 Sc only. 

fer that m 6 conc ® ntration is known, the catalytic heat trans- 

followsf transfer due to diffusion, q d ^ w , is obtained as 


^■d,w Pwkw c i,w^h - 


kwpMg \ 

. RT W ) 


Ah° 


m w c i,e \ 

l+mvci, e y 


( 68 ) 


the total heat transfer is given by 


q w " ^d,w + ^c,w (69) 

where q c/w - is the heat transfer due to ordinary convection, values of 
which can be obtained elsewhere, such as in references 7 and 9. 


Axisymmetric Conical Body With Spherical Nose 


In this section, equation (52) will be first solved numerically 
for the axisymmetric conical body with spherical nose shown in fig- 
ure 1 (b) . An approximate solution will also be obtained by the use of 
the local similarity concept and the two solutions will be compared. 
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For simplicity, the following approximations used in reference 15 
are used here for the inviscid flow properties. 

X = cos2X , u e = u^X for xSe (7< 

Po 

~ = cos 2 e , u e = uJ3 o e for X = 0 (7- 

Po 

The streanrwise variable s becomes for the bcdy 

S = Cu^PqL 3 f X Sin2x co S 2X dX for X = 0 (7 

J o 


S = CuJ3 0 L 3 ^ / x sin 2 X cos 2 X dX 


+ 0 CO£ 


S 20 T X [sin 0 + (x- 0 )cos 0] 2 dx| for X ^ 0 U 


also the function 0 becomes 


X sin 2 X cos X dX 


p Q X sin^X 


for X = e 


X sin 2 X cos 2 X dX + 


n X .,2 

0 cos 2 0 / [sin 0 + (X-0 )cos 0] dX 

Je 


P o 0 2 [sin 0 + (X-0)cos 0]' 


for X ^ © (75) 


Now to remove a flight condition from the shape factor let 

5 » VPo ® 


and also define 


ScG 

^ B i Bl 


(77) 
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Now with the aid of equations (70) through (77), equation (52) becomes 

x\ 

/ l+*%c i L e 


cbv . 1 - "» - l. 1 + ^ x + x * 


sin X cos^ 


1 + 


Kj® 


(l+m w c 1 ^ e ) ; 


$ 2 X 

cos 2 X 


for X ^ 6 ( 78 ) 


dm 

dX 


■w 


1 - 1% - 2K;]$ -l 

*i + . , 1 ® 2q 

L m w 

L A B 2 Sc cos 0L sin 0 + (X-0)cos 0) J 

l+myc 1 e 

2 

[1 + *3* - 

O 2 0 

1 L- 

(l+rnyC!^) 2 ] 

cos 2 0 


for X = e ( 79 ) 

The solution of equation ( 78 ) for X = 0 (stagnation point) can be 
obtained from the general local similarity solution of equation ( 35 ) 

This solution is for T = 0 u 


m (2 + 1 c i,e) + J (o, + 1 - c 1)6 ) 2 + 4 c]^ e 


where 

ft - Kjl ( 1 + — ^ 0 
V A'B 2 Sc/ 

The value of ft at the stagnation point is found by noting 

lim I = l /2 
X-O 

Equation ( 80 ) was found to yield stagnation point values which check 
within 5 percent of the exact solution given in reference 6 . 

For the region downstream of the stagnation point equations ( 78 ) 
and (79) were integrated numerically for Sc = 0.72 and several values 
of G and c^e, and with 9 equal to 50° and 70°. These results are 
shown m figure 4. The solutions for G = 0 and G = 00 represent the 
two extremes of catalytic activity. 


(80) 

(81) 
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Now let us examine the possibility of using the local similarity 
solution everywhere along the body. In order no estimate the accuracy 
of the local similarity solution, the hypersonic pressure and equilibn 
heat- transfer distributions (obtained from ref. 15) are shovn below for 
tfte body in question. 



Sketch (a) 

Since the curves ore quite alike except along the conical afterbody for 
I” 50°, equation ( 80 ) is expected to be generally quite satisfactory 
in the determination of the atom wall concentration. 

From figure 4 we see that the variation of m w with X becomes 
less pronounced as G approaches either of the two extremes. The error 

introduced by a local similarity assumption should t*®*®^®*® ere^it 
at intermediate values of G, and if its use can be justified there, it 

should be valid for all values of G. 

Figure 5 shows the catalytic heat transfer based on the numerical 
solution and the local similarity solution for an intermediate value of 
the reaction parameter G = 5-0 and c x , e = 1-0. The agreement between 
the two results is strikingly good, with some deviation along the conical 
afterbody 6 for 6 = 50° as expected. It seems that the local similarity 
solution^ is quite adequate for the calculation of catalytic heat transfer 

to this body. 
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APPLICATION TO HYPERSONIC FLIGHT 
Altitude Range of Application 


As flight altitude is increased, the boundary layer surrounding a 
hypersonic vehicle passes through a regime in which homogeneous recom- 
bination is significant into one where this reaction is frozen and, 
finally, with very low density, into a regime where the classical 
boundary- layer theory is no longer valid. The present analysis is 
applicable in the intermediate regime, that of near frozen boundary- 
layer flow, where catalytic recombination is probably most significant 
in the calculation of heat transfer. 

Heat transfer to blunt bodies with nonequilibrium boundary layers 
was considered in reference 17 > where it is shown that the boundary 
layer at the stagnation point becomes practically frozen at about 

200.000 feet altitude, depending, of course, on the nose radius and wall 
temperature. Also, for a particular nose radius and wall temperature, 
the boundary layer on the conical afterbody will be frozen at consider- 
ably lower altitudes. 

As the Reynolds number decreases with increasing altitude, the vor- 
ticity in the inviscid flow generated by the bow shock becomes compar- 
able to the vorticity in the boundary layer. This is called the vorti- 
city interaction regime, in which it is necessary to modify the classical 
boundary- layer theory. 

Reference l8 indicates that the boundary-layer concept is valid near 
the stagnation point of a blunt body whenever M^/Re^ « p^/p 0 e (b/A) . 
Consider a hypersonic body of 1 foot radius flying at an altitude of 

250.000 feet. The ratio, M^/Re^, is about one order of magnitude smaller 
than Poo/Po,ej which suggests that boundary- layer theory can still be 
approximately applied at this altitude. The effect of vorticity inter- 
action on ordinary convective heat transfer to the stagnation point was 
analyzed in reference l8 for a compressible inert gas. It can be cal- 
culated from this analysis that the vorticity interaction effect on heat 
transfer to the stagnation point of a sphere at the altitude of 250,000 
feet is to increase the heat transfer about 10 percent. Since the vor- 
ticity in the boundary layer increases rapidly as the gas flows around 
the body, the altitude at which the vorticity interaction effect becomes 
important on the conical afterbody should be considerably higher. The 
present analysis should therefore be applicable between about 200,000 
feet and 250,000 feet, a range which may be extended considerably in both 
directions if one is primarily interested in the conical afterbody. 
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Discussion of Gas Mode] 


Air, for the present flight range, is more validly described as a 
mixture of oxygen and nitrogen than as a single diatomic gas. It is felt 
however, that since the transport properties are quite alike, any impor- 
tant difference, as regards catalytic heat trar.sfer, would be in the 
reaction rates and heats of reaction. 

The conservation of species equation when written independently for 
oxygen and nitrogen would be identical and wher. combined, would give 
equation (35)* Now, if k v were the same for both oxygen and nitrogen, 
the combined catalytic boundary condition woulc. be the same as that for 
a single diatomic gas, equation (48). The assumption of equal values of 
k w seems reasonable in view of the limited knc wledge of heterogeneous 
reaction rates (ref. 6), and is made here. 

With these limitations, equal transport properties and equal values 
of k v , the combined atom wall concentration, nay be calculated using the 
present method. 

If the solution had been carried out independently for oxygen and 
nitrogen, with the above assumptions, the reduced concentration profile 
would be the same for both since the equations are identical with identi- 
cal boundary conditions. The heat of recombination at the wall Ah 0 may 
therefore, in the present analysis, be obtainec. as a weighted average 
based on the atom concentrations at the boundary- layer edge. 


The Effect of Catalytic Activity on Hypersonic Heat Transfer 


For a given hypersonic blunt body, the catalytic heat transfer is 
largely controlled by the reaction parameter G as can be seen from the 
general equation (52) . Approximate values of G for two typical surface 
materials are shown in figure 6 for the range of flight conditions of 
interest. The calculations are for the blunt body shown in figure 1 with 
a nose radius of 1 foot. The values of k w are taken from reference 6. 

The atom diffusion rate into the wall is approximately proportional 
to the driving potential, e - c x v . Therefore, for m w = 0, q^ v is 
maximum and the total heat transfer would be about the same as that for 
an equilibrium boundary layer. The proportion of dissociation energy 
transferred to the wall will decrease with ( 1— r i w ) . Now let us consider 
a metallic oxide surface, which is quite highly catalytic. It is seen 
from figures 4 and 6 that at lower altitudes ('.ess than 150,000 ft) the 
total heat transfer, q v , will be about as largo as equilibrium heat trans- 
fer. At the altitude of 250,000 feet, however, the catalytic heat trans- 
fer will be reduced by about 50 percent. Since, at high speeds, as much 
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as 75 percent of the total energy could be associated with dissociation 
the total heat transfer would be reduced to about 65 percent of equilib- 

h °f Tt' SVen f ° r thlS hlgh1 ^ “talytlc surface. My surface 
material with a lower catalytic efficiency, such as pyrex, would, of 

course, reduce the heat transfer to a much lower value. 


CONCLUDING REMARKS 
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axi^^+ Cata i y ^ iC heat transfer to the surfaces of two-dimensional and 
axisymmetnc bodies with finite catalytic activities was analyzed. 

An integral method was used to obtain catalytic heat-transfer pre- 
dictions to (1) a flat plate in frozen compressible flow with and with- 

nosP T 1 ^ ^ ^ axis y mmetric conical body with a spherical 

1 frozen hypersonic compressible flow without transpiration. 

The solution for a flat plate was derived in closed form. 

aho ^ n that ^ a body whose pressure and equilibrium heat 
transfer distribution curves are similar, the catalytic heat transfer 
to the body can be calculated quite accurately by the use of a simple 
oca similarity concept. For the axisymmetric conical body with spher- 
ical nose, the pressure and the equilibrium heat-transfer distribution 
patterns were found to be quite similar. A comparison of the catalytic 

tha^oMaSd r ?I 10n baS6d ° n the Simple l0Cal simila rity concept with 
th t obtained by the numerical solutions showed that the maximum error 

was in the order of 5 percent for this body. 

Pt was shown that at high flight altitudes, the dissociation energy 
would not be recovered completely even with highly catalytic surfaces, 
such as a metallic oxide. The total heat transfer, therefore, may be 

considerably less than the equilibrium value for most engineering 
incite 0 r* x £lX s * 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 22, i960 



26 


references 


1 Chambre Paul L., and Acrivos, Andreas: On Chemical Surface Reactions 

in Laminar Boundary Layer Flows. Jour. Appl. Phys., vol. 27, 
no. 11, Nov. 1956, PP. 1322-1328. 

. - nrir > Phambre Paul L.: Laminar Boundary Layer Flows 

2. Acrivos, Andreas, and Chambre, Paul n . nQ< 6 June 1957, 

with Surface Reactions. Ind. Eng. onem., y> 

pp. 1025-1029. 

, * _ , T . nn Phpmical Surface Reactions in Hydrodynamic 

3 - Tpl Sci R-T-. *1 vol. 6, no. 2-3, 1956, pp. 97-113- 

I _ -^1 T? . Phpmicallv Frozen Boundary Layers with Catalytic 

4 - ^°a,rf ace^eact ion . ‘SSTSo/U. *i., Si. *, - 5, May 1959, 

PP* 281-286. 

a Rosner Daniel E.: Diffusion and Chemical Surface Catalysis in Flou 

Systems. Tech. Pub. Ho. 14. AeroChem. Research Labs. Inc., 
Princeton, N. J., Sept. 1959 • 

6. Goulard, Robert J.: On Catalytic Recombination Rates in 

Stagnation Heat Transfer. Jet Propulsion, vol. 28, no. 11, 

1958, pp. 737-745- 

7 Chung, Paul M.: Shielding Stagnation Surfaces of finite Catalyt 10 

^Activity by Air Injection in Hypersonic Flight. NASA TN ^27, 

1959. 

8 . Ch^g, Paul M., and person 

pSe! le AS°Lu f r.: vol. 30, ro. 3, Mar. i960, PP . 262-264. 

_ . » TT. a+ Tranqfer With Mass Addition And Chem- 

9 * Lee ical L ReaSions! nV ?hird AGARD Combustion and Propulsion Panel Collo- 
JuS, Palermo, Sicily, Mar. 17-21, 1958. Pergamon Press, N. Y., 
pp. 451-498. 

10 ' Ke TrLsfer°Around^Blunt 41 BoSes*in a Dlssociated^lr.^*JourT I1 Aero /Space 
Sci., vol. 26, no. 7, July 1959, PP- >> 2 -. -430. 

11. Rosner, Daniel E.: Boundary Conditions F>r the Flov of Multicomponent 

Gas. Jet Propulsion, vol. 28, no. 8, pc. 1, Aug. 1950, PP- 757 70 

12. Sc ala, Sinclaire M : 

Surfaces Having Finite Catalytic Eilici y 
gress of Appl. Mech. 1958, PP- 799-806. 


A 

4 
2 

5 



< -=t oj m 


27 


13. Cutting, John C. , Fay, James A., Hogan, William T., and Moffatt, W. 
Craig: A Non-catalytic Surface for Dissociated Combustion Gases. 

Preprint, Confer, on Physical Chem. in Aero, and Space Flight, 

Sept. 1-3, 1959. 

lb. Morduchow, Morris: Analysis and Calculation by Integral Methods of 

Laminar Compressible Boundary Layer With Heat Transfer and With and 
Without Pressure Gradient. NACA Rep. 1245, 1955. 

15. Lees, Lester: Laminar Heat Transfer Over Blunt-Nosed Bodies at 

Hypersonic Flight Speeds. Jet Propulsion, vol. 26, no. 4, Apr. 
1956, pp. 259-269. 

16. Low, George M: The Compressible Laminar Boundary Layer with Fluid 

Injection. NACA TN 34o4, 1955. ' 

17 • Chung, Paul M., and Anderson, Aemer D.: Heat Transfer Around Blunt 

Bodies with Nonequilibrium Boundary Layers. Proc. i960 Heat Trans- 
fer and Fluid Mechanics Institute. 

l8. Hayes, Wallace D. , and Probstein, Ronald F.: Hypersonic Flow Theory. 

Academic Press, New York, 1959, pp. 370-386. 



28 


ui ro 4=- > 


























